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Abstract. Direct and indirect dark matter detection relies on the scattering of tiie dark matter 
candidate on nucleons or nuclei. Here, attention is focused on dark matter candidates (neutralinos) 
predicted in the minimal supersymmetric standard model and its constrained version with universal 
input soft supersymmetry-breaking masses. Current expectations for elastic scattering cross sections 
for neutralinos on protons are discussed with particular attention to satisfying all current accelerator 
constraints as well as insuring a sufficient cosmological relic density to account for the dark matter 
in the Universe. 



1. Introduction 

^ . Minimal supersymmetric theories with R-parity conservation are particularly at- 

tractive for the study of dark matter as they predict the existence of a new stable 
particle which is the lightest R-odd state (the LSP). Furthermore, for parameter 
values of interest to resolve the gauge hierarchy problem, the LSP has an annihi- 
I lation cross section which yields a relic density of cosmological interest (Ellis et 

al., 1984). Recent accelerator constraints have made a great impact on the available 
parameter space in the MSSM and in particular the constrained version or CMSSM 
in which all scalar masses are assumed to be unified at a grand unified scale (Ellis 
et al., 1998; EUis et al., 2000a; EUis et al., 2001b; EUis, Nanopoulos, & Olive, 
2001). In addition, significant progress has been made concerning the relic density 
calculations (Ellis et al., 2000b; Ellis et al., 2001a; Gomez, Lazarides, & Pallis, 
2000). These issues were discussed in detail in John Ellis' contribution and will 
only be touched on briefly here. 

The main impact of the null searches, particularly at LEP, is in the increase in the 
lower limit to the LSP mass as well as the rest of the sparticle spectrum. This has the 
unfortunate effect of lowering the elastic scattering cross sections for neutralinos 
on protons, making direct detection experiments more difficult. Nevertheless, we 
have now entered a period where accelerator constraints will be at a lull due to the 
transition from LEP to the LHC and the time required for run II at the Tevatron to 
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acquire the needed luminosity*. Therefore, during the next several years we can be 
hopeful that direct detection experiments will improve to make meaningful inroads 
to the supersymmetric dark matter parameter space. 

In this contribution, we will discuss the current status of the expected neutralino- 
proton elastic cross sections in the MSSM and CMSSM. These results will be 
applied to the possibility of direct dark matter detection. 

As noted earlier, we will restrict our attention to regions of parameter space for 
which the relic density of neutralinos is of cosmological interest. For an age of the 
Universe t > 12 Gyr, there is a firm upper bound on the relic density fi^/i^ < 0.3, 
where is the fraction of critical density in the form of neutralinos, %, and h 
is the Hubble parameter in units of 100 km/s/Mpc. This limit represents a strict 
cosmological bound on the supersymmetric parameter space. We also focus our 
discussion on the parameter values which lead to relic densities with fi^/i^ > 0.1. 
While this does not place any bound on supersymmetry, it is a reasonable require- 
ment for dark matter candidates. Neutralinos with a lower density could not be the 
dominant form of dark matter in our Galaxy, and therefore detection rates would 
necessarily be suppressed. 



2. The MSSM vs. The CMSSM 

As discussed here by John Ellis, the neutralino LSP is the lowest-mass eigenstate 
combination of the Bino B, Wino W and Higgsinos Hi^2, whose mass matrix A'^ 
is diagonalized by a matrix Z: diag{m^^^ ) = Z*NZ~^. The composition of the 
lightest neutralino may be written as 

X = Z^iB + Z^2W + Z^sHi + Z^^H2 (1) 

We assume universality at the supersymmetric GUT scale for the gauge couplings 
as well as gaugino masses: Mi^2,3 = 'mi/2, so that Mi ~ | tan^ 9wM2 at the 
electroweak scale (note that this relation is not exact when two-loop running of 
gauge sector is included, as done here). 

We also assume GUT-scale universality for the soft supersymmetry-breaking 
scalar masses mo of the squarks and sleptons. In the case of the CMSSM, the 
universality is extended to the soft masses of the Higgs bosons as well. We fur- 
ther assume GUT-scale universality for the soft supersymmetry-breaking trilinear 
terms ^o- In the MSSM, we treat as free parameters mi/2 (we actually use M2 
which is equal to mi/2 at the unification scale), the soft supersymmetry-breaking 
scalar mass scale mo (which in the present context refers only to the universal 
sfermion masses at the unification scale), Aq and tan /3. In addition, we treat /x 
and the pseudoscalar Higgs mass tua as independent parameters, and thus the two 

* We can however, expect improvements in tire uncertainties in rare B decays and tlie measure- 
ment of the anomalous magnetic moment of the muon, both of which will have an impact on the 
allowed supersymmetric parameter space. 
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Higgs soft masses mi and m2, are specified by tlie electroweak vacuum conditions, 
which we calculate using mt = 175 GeV. In contrast, in the CMSSM, mi and 
7712 are set equal to mo at the GUT scale and hence fi (up to a sign) and m^ 
are calculated quantities, their values being fixed by the electroweak symmetry 
breaking conditions. 



3. Elastic Scattering Cross Sections 



The MSSM Lagrangian leads to the following low-energy effective four-fermion 
Lagrangian suitable for describing elastic x-nucleon scattering (Falk, Ferstl, & 
Olive, 1999): 



(2) 



This Lagrangian is to be summed over the quark generations, and the subscript 

i labels up-type quarks {i = 1) and down-type quarks (i = 2). The terms with 
coefficients an, a4i,a5i and ctei make contributions to the elastic scattering cross 
section that are velocity-dependent, and may be neglected for our purposes. In fact, 
if the CP-violating phases are absent as assumed here, 05 = ag = (Falk, Ferstl, 
& Olive, 2000; Chattopadhyay, Ibrahim & Nath, 2000). The coefficients relevant 
for our discussion are the spin-dependent coefficient, 02, 
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and the spin-independent or scalar coefficient, 0:3, 
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Yi = ryn [-g Z^, + gTs.Z^^) + ^^^^^ 



TXT — * ^"^9i-^x5-i * I ry* 

Vi = r?2*2^?^ + %*i i^g'Z^, + gn,Z^2 ) (5) 



2mwBi 

where yi, T^i denote hypercharge and isospin, and 

^li = Zy.3{Z^4) , 52i = Z^4{-Zxz), 

Bi = svn [3 {cos jS) , = cos /3(— sin /3), 

Cj = sina(cosa) , = cosa(— sina) (6) 

for up (down) type quarks. We denote by < "m-Hi the two scalar Higgs 

masses, and a denotes the Higgs mixing angle. Finally, the sfermion mass-squared 
matrix is diagonalized by a matrix rj: diag{m\,m^) = r)M'^r]~^, which can be 
parameterized for each flavour / by an angle 9f: 

cos 9f sin 0/ ^ ^ f r/ii rju , 



-SmUf cosy/ J y 7/21 r)22 

The spin-dependent part of the elastic x-nucleus cross section can be written as 

^2 = — G|.m2AV( J + 1) (8) 

TT 

where rrir is again the reduced neutraUno mass, J is the spin of the nucleus, and 

A = ^{ap{Sp) + an{Sn)) (9) 
where 

The factors A^-''''^'' parameterize the quark spin content of the nucleon. A recent 
global analysis of QCD sum rules for the gi structure functions (Mallot, 1999), 
including C(a^) corrections, corresponds formally to the values 



A^P) = 0.78 ± 0.02, aJ'^ = -0.48 ± 0.02 

A^) = -0.15 ± 0.02 (11) 
The scalar part of the cross section can be written as 

a, = ^[Zfp + {A-Z)Uf (12) 

TT 
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where 



rrir, ^-^ ^ rrig 27 ^ rria 

q=u,d,s c,b,t 

and /„ has a similar expression. The parameters f^^ are defined by 

"^pfrq = (plmqqqlp) = mqBq (14) 

whilst f!^l = 1 - Y.q=u,d,s frq (Shifman, Vainshtein, & Zakharov, 1978). We 
observe that only the products niqEq, the ratios of the quark masses ruq and the 
ratios of the scalar matrix elements Bq are invariant under renormalization and 
hence physical quantities. 

We take the ratios of the quark masses from (Leutwyler, 1996): 

77? TTt 

— = 0.553 ± 0.043, — = 18.9 ± 0.8 (15) 

md nid 

In order to determine the ratios of the Bq and the products mqBq we use infor- 
mation from chiral symmetry applied to baryons. Following (Cheng, 1989), we 
have: 

By, - Bs m^o + ms- -mp-mn 

Bd - Bs ms+ +m,Y,- -mp- mn 

Substituting the experimental values of these baryon masses, we find 

z = 1.49 (17) 

with an experimental error that is negUgible compared with others discussed below. 
Defining 

Bd + Bu 



we then have 

Bd 2 + ((z - 1) X y) 



S„ 2 X z - ((z - 1) X y) 



(19) 



The experimental value of the vr-nucleon a term is (Gasser, Leutwyler & Sanio, 
1991; Knecht, 1999): 

a = ^(m„ + md) X + = 45 ±8 MeV (20) 

and octet baryon mass differences may be used to estimate that (Gasser, Leutwyler 
& Sanio, 1991; Knecht, 1999) 

a = , ^° , : (70 = 36 ± 7 MeV (21) 

(1 -y) 
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The larger value of cr = 65 MeV (Olsson, 2000; Pavan et al., 1999) considered 
by (Amowitt, Dutta, & Santoso, 2000) leads to scattering cross section which are 
larger by a factor of about 3. Comparing (20) and (21), we find a central value of 
y = 0.2, to which we assign an error ±0.1, yielding 

— = 0.73 ± 0.02 (22) 

Bu 

The formal error in y derived from (20) and (21) is actually ±0.2, which would 
double the error in B^/Bu. We have chosen the smaller uncertainty because we 
consider a value of y in excess of 30% rather unlikely. 

The numerical magnitudes of the individual renormalization-invariant products 

(p) 

TUqBq and hence the may now be determined: 

f^l = 0.020 ± 0.004, /^^J = 0.026 ± 0.005 

f!^} = 0.118 ± 0.062 (23) 

where essentially all the error in /^^ arises from the uncertainty in y. The corre- 
sponding values for the neutron are 

f^^ = 0.014 ± 0.003, f^J = 0.036 ± 0.008 

= 0.118 ±0.062 (24) 

It is clear already that the difference between the scalar parts of the cross sections 
for scattering off protons and neutrons must be rather small. 



4. Results 

We begin by discussing the results for the CMSSM (Ellis, Ferstl, & Olive, 2000). 
For fixed tan f5 and sign of ^, we scan over experimentally and cosmologically 
allowed regions in the mi/2 — friQ plane. Results here are shown for Aq = 0. The 
combination of the cosmological constraint Q.h? < 0.3 and the constraint from the 
Higgs mass, mj^j > 113 GeV, eliminates low values of tan/3 < 5 (Ellis et al., 
2001b). For the value of tan/3 = 10, we show in Figure 1 the elastic scattering 
cross section for spin-dependent (a,b) and scalar (c,d) processes as a function of 
the neutralino mass. Although it is barely discernible, the thicknesses of the central 
curves in the panels show the ranges in the cross section for fixed that are 
induced by varying mo. At large m^^ where coannihilations are important, the range 
in the allowed values of mo is small and particularly little variation in the cross 

section is expected. The shaded regions show the effects of the uncertainties in 

(v) (v) 
the input values of the A^- (11) (a,b) and in the frf' (23) (c,d). For the results of 

analogous analyses, see (Accomando et al, 2000; Corsetti & Nath, 2000; Bottino 

et al., 2001). 
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In addition we show the constraint coming from mj^j > 113 GeV for = 
which restricts one to relatively large neutralino masses. For the cases where > 
0, there is a potential upper limit to the neutralino mass coming from the recent 
BNL E821 experiment ( Brown et al, 2001), which reports a new value for the 
anomalous magnetic moment of the muon: — 2 = 2 x that is apparently 
discrepant with the Standard Model prediction at the level of 2.6 a. This limit is 
also displayed on Figure 1 (ElUs, Nanopoulos, & Olive, 2001; Amowitt, Dutta, & 
Santoso, 2001). For tan /S = 10, fi > 0, the theory is quite predictive in both the 
LSP mass and scattering cross section. No value of tan /3 is compatible with the 
BNLE821 result if |L( < 0. 

The scalar cross section is, in general, more sensitive to the sign of /x than is the 
spin-dependent cross section. Notice that, in Figure Ic for tan (3 = W and /U < 0, 
there is a cancellation. Higgs exchange is dominant in 03 and for < 0, both 
Zy^3 and are negative, as is the Higgs mixing angle a. Inserting the definitions 
of 5ii(2t), we see that there is a potential cancellation of the Higgs contribution to 
as for both up-type and down-type quarks. Whilst there is such a cancellation for 
the down-type terms, which change from positive to negative as one increases m-^, 
such a cancellation does not occur for the up-type terms, which remain negative 
in the region of parameters we consider. The cancellation that is apparent in the 
figure is due to the cancellation in 03 between the up-type contribution (which is 
negative) and the down-type contribution, which is initially positive but decreasing, 
eventually becoming negative as we increase m^. 

At higher values of tan one can in principle expect larger elastic cross sec- 
tions (Accomando et al., 2000; Lahanas, Nanopoulos, & Spanos, 2000). In Figure 
2, we show the spin dependent (a,b) and the scalar (c,d) cross section for tan (3 = 
35, fi < (a,c) and for tan (3 = 50, fi > (b,d). In this figure, lower values 
of niy^ have been cut off (and are not shown) due to the constraint imposed by 
measurements ofb^sj. 

As was discussed in detail in (Ellis et al., 2001a), a new feature in the mi/2 — rno 
plane with acceptable rehc density appears at large tan/?. At large mi/2 ^ 1000 
GeV, it becomes possible for neutraUnos to annihilate through s-channel Hi or 
pseudoscalar. A, exchange. In fact there is a slice in the plane where 2m^ 'mHi,A 
and the relic density becomes uninterestingly small. At smaller and larger 1111/2 
surrounding this pole region, there are regions where the relic density falls in the 
desired range. This leads to two separate regions in Figure 2 at lower m^^. The third 
region in Figure 2 at higher m^^ corresponds to the cosmological region allowed by 
coannihilation (Ellis et al., 2000b). For more further details on i^i, A-pole and 
coannihilation, see the contribution of John Ellis in these proceedings. As in the 
case tan /? = 10, < 0, the scalar cross section at higher tan (3 also exhibits the 
cancellation feature discussed above. However, because the cosmological regions 
are multivalued in mo as a function of m 1/2 > the cancellation occurs at a different 
value of for each of three regions just discussed. This leads (unfortunately) to 
a broad region in the cr — m^^ plane where the cross section is very small. 
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Figure 1. (a,b): The spin-dependent cross section for the elastic scattering of neutralinos on protons 
as a function of the LSP mass for tan/3 = 10. The central curves are based on the inputs (11), 
and their thicknesses are related to the spreads in the allowed values of mo. The shaded regions 
correspond to the uncertainties in the hadronic inputs (11). (c,d): The spin-independent scalar 
cross section for the elastic scattering of neutralinos on protons as a function of the LSP mass for 
tan f3 = 10. The central curves are based on the inputs (23), their thicknesses are again related to 
the spread in the allowed values of mo, and the shaded regions now correspond to the uncertainties 
in the hadronic inputs (23). The supplementary lower limits imposed on in this and the next 
figure reflect improvements in the LEP lower limit on mu, and the upper limits for p > are due to 
Qfj, — 2, which is incompatible with p < 0. 
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Figure 2. As in Figure 1 for tan 13 = 35. Here the three distinct regions correspond to the two sides 
of the Hi, A annihilation poles, and the coannihilation region at higher values ofm^. 

In the MSSM, in addition to scanning over the gaugino and sfermion masses 
at fixed tan (i, one can treat ^ and the pseudoscalar mass niA as free parameters 
as well. In (Ellis, FerstI, & Olive, 2001), we performed a scan over the following 
parameter space: < mo < 1000; 80 < < 2000; 80 < Ms < 1000; < 
mA < 1000; -1000 < yl < 1000. Of the 90,000 (70,000) points scanned for 
tan /3 = 10 and /i > (/i < 0), only 6208 (4772) survived all of the experimental 
and cosmological constraints. In the CMSSM, the LSP is nearly always predicted 
to be Bino of very high purity. However, in the MSSM, when |/i| < M2, the LSP 
may have a dominant Higgsino component. In these cases, coannihilation (Griest 
& Seckel, 1991) greatly suppresses their relic density and when combined with 
the experimental constraints on the parameter space, Higgsino dark matter can be 
excluded as a viable option (Ellis et al., 1998; Ellis et al., 2000a). 

The LEP chargino and Higgs cuts remove many points with low and/or 
large elastic scattering cross sections. The sfermion mass cut is less important. The 
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constraint that x be the LSP removes quite a large number of points, populated 
more or less evenly in the cross section plots. There is a somewhat sparse set of 
points with very small cross sections which give some measure of how low the 
cross section may fall in some special cases. These reflect instances where particu- 
lar cancellations take place, as discussed above. The lower boundary of the densely 
occupied regions offers an answer to the question how low the elastic scattering 
cross sections may reasonably fall, roughly a ~ 10~^ pb for the spin-dependent 
cross section and ~ 10~^° pb for the spin-independent cross section. 

Our resulting predictions for the spin-dependent elastic neutralino-proton cross 
section for tan (3 = 10 are shown in Figure 3(a,b), where a comparison with the 
CMSSM is also made. The raggedness of the upper and lower boundaries of the 
dark (blue) shaded allowed region reflect the coarseness of our parameter scan, 
and the relatively low density of parameter choices that yield cross sections close 
to these boundaries, it should be noted that the low values of m^^ in these plots, that 
yield relatively high spin-dependent cross section, have now been excluded by im- 
provements in the Higgs mass hmit. As increases, the maximum allowed value 
of (J spin decreases, though not as rapidly as in the previous CMSSM case (ElUs, 
Ferstl, & Olive, 2000). The hadronic uncertainties are basically negUgible for this 
spin-dependent cross section, as seen from the light (yellow) shading. 

The analogous results for the spin-independent elastic neutralino-proton cross 
section are shown in Figure 3(c,d), where comparisons with the CMSSM case are 
also made. We see a pattern that is similar to the spin-dependent case. For small 
m^, the spin-independent scalar cross section, shown by the dark (blue) shaded 
region, may be somewhat higher than in the CMSSM case, shown by the (red and 
turquoise) diagonal strip, whilst it could be much smaller. For large m^, the cross 
section may be rather larger than in the CMSSM case, but it is always far below the 
present experimental sensitivity. Overall, we note that the hadronic uncertainties, 
denoted by the light (yellow) bands, are somewhat larger in the spin-independent 
case than in the spin-dependent case. 



5. Conclusions 

As one can see from scanning the figures, the predicted elastic scattering cross 
section in the CMSSM and in the more general MSSM, are relatively small. For 
the spin-dependent processes, the cross sections fall in the range a ~ 10~^ — 10~^ 
pb, whereas for the scalar cross sections, we find a < 10~^ pb with an uncertain 
lower limit due to possible cancellations. These should be compared with current 
sensitivities of existing and future experiments (Gaitskell & Mandic, 2001). The 
UKDMC detector is sensitive to a > 0.5 pb for the spin-dependent cross sec- 
tion. DAMA and CDMS are sensitive to o" > 2 x 10~^ pb for the scalar cross 
section. This is close to the upper limits we find for reasonable supersymmetric 
models. The future looks significantly brighter. When CDMS is moved to the 
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Figure 3. As in Figure 1. The main (blue) shaded regions summarize the envelopes of possible values 
found in our scan, for points respecting the LEP constraints, discarding points with Q,^h? > 0.3, 
and reseating points with Q.^h? < 0.1. The small light (yellow) shaded extensions of this region 
reflect the hadronic matrix element uncertainties. The concave (red and turquoise) strips are those 
found previously assuming universal Higgs scalar masses (Ellis, Ferstl, & Olive, 2000j. 
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Soudan mine, its sensitivity will drop to between 10~^ and 10"'' pb and GENIUS 
claims to be able to reach 10~^pb. At those levels, direct detection experiments 
will either discover supersymmetric dark matter or impose serious constraints on 
supersymmetric models. 
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